Segregation of Polymers in Confined Spaces 
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We investigate the motion of two overlapping polymers confined in a 2d box. A statistical model is 
constructed using blob free-energy arguments. We find spontaneous segregation under the condition: 
L > _R||, and mixing under L < Tiy, where L is the length of the box, and _R|| the polymer extension 
in an infinite slit. Segregation time scales are determined by solving a mean first-passage time 
problem, and by performing Monte Carlo simulations. Predictions of the two methods show good 
agreement. Our results may elucidate a driving force for chromosomes segregation in bacteria. 



Biopolymers have evolved to function in crowded and 
confined environments [ij. Under these conditions, ex- 
cluded volume effects and geometrical confinement com- 
pete with entropy to yield unique structures and dynam- 
ical processes [1, d, 13, @|- Examples include the trans- 
portation of proteins through a membrane channel and 
endoplasmic reticulum, and the replication and segrega- 
tion of highly compacted chromosomes during cell divi- 
sion [i,0,li]- 

The mechanism underlying the process of chromo- 
somes segregation in bacterial cells is still unclear 0, 
H, [l3| ■ Recent results of molecular dynamic simulations 
indicate that entropic driving forces, in the absence of 
motor proteins, can describe some important features 
of chromosome segregation in C.crescentus and E.coli 
0) El • Analytic calculations of the dynamics in an open 
cylinder predicts a constant segregation velocity [13] . Ex- 
perimental studies of the sequential movement of chromo- 
somal loci during replication, however show that the seg- 
regation process is quick and inhomoeeneou s [ll . IsL 1^ . 
Analysis of segregation in a closed geometry [15[ indicates 
that the shape of the confining box determines whether 
polymers segregate or remain mixed. A natural question 
that arises is whether the inhomogeneous dynamics ob- 
served in experiments is a consequence of the shape of 
the cell. 

In this letter, we analyze the motion of two identical, 
self-avoiding chains, confined in a 2d rectangular box. 
Combining scaling theory, stochastic models, and Monte 
Carlo simulations, we show the existence of a transition 
from segregation to mixing with change of aspect ratio. 
The dynamics of segregation is shown to be inhomoge- 
neous, with two distinct regimes. The segregation time 
depends sensitively on the geometry, and interestingly, 
exhibits a minimum at a specific aspect ratio that de- 
pends on the density. 

Two identical self-avoiding chains with N monomers 
are confined to a 2d box of width W and length L, as il- 
lustrated in Fig[T]Ja). The width is chosen to be less than 
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FIG. 1: (Color online) (a) Blob structure of two polymers con- 
fined to a rectangular box. _R|| is the extension of confined 
polymer (cf. text). The illustration is schematic and only 
the size of the blob, ^uot, is relevant. (6) Illustration of box 
partitioning during segregation. Region I is overlapping and 
region II, III are single-chain. 



blob size, ^uob, is determined by the competition between 
confinement, excluded- volume effects, and entropy. For 
length scales small compared to S,biobi confinement effects 
are screened and self-avoidance effects dominate, but for 
length scales larger than the blob size, polymer confor- 
mations are strongly affected by the confining geometry. 

For a single chain confined in a rectangular box, assum- 
ing homogeneous distribution of monomers 1^, 
leads to 
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(1) 



Rg, the radius of gyration of the free polymer. A scaling 
theory of the free energy of confined polymers is based 
on the blob pictures with each blob contributingfcsT to 
the free energy of a polymer (Fig[Ha)) [3, E3|- The 



where = | is the Flory exponent in 2D, and the length 
is measured in units of the Kuhn length [l^, [l3|. The 
width of the box provides an upper bound for ^uob, which 
defines a length scale R\\ cx NW~^ 16] that measures the 
extension of the polymer in the longitudinal direction. 
In a rectangular box, if L » blobs occupy only a 
small portion of the confining space, and the situation 
is similar to an open slit. In the opposite limit, two- 
dimensional confinement should significantly affect the 
properties of the chain. We introduce a dimensionless 
parameter A = — r — , to measure the infiuence of the 
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closed geometry; A > 1 if L < i?|| , and A < if X > 2i?||. 
Since each blob contributes ~ ksT = 4 to the free energy 



of a single polymer [B, [l3l, F{L, W, N) 



pF{L, W, N) oc Nuob 



if A > 1 
if A < 1 



(2) 



where Nbiob is the number of blobs in a single confined 
polymer. 

In order to model two overlapping polymers, as illus- 
trated in Figllja), the blob picture can be extended, as- 
suming (a) that the density is uniform, and (b) that 
is the same as that of a single polymer (limit of infinites- 
imal thickness of chains). The monomer linear density 
for two overlapping polymers is: 



^ ifA>l 
if if^<l 



(3) 



During the process of segregation, the box space is natu- 
rally partitioned into three regions: Region I of length I 
where the two polymers overlap, and Regions II and III 
that are occupied by single-polymer segments (Fig[T]^b)). 
By symmetry, regions II and III are equivalent. The total 
free energy of the polymers is a function of the overlap 
distance l,0<l<Rf F{1) = Fo + 2Fs where Fo, F, de- 
note the free energy in the overlapping and single-chain 
regions, respectively. If the dominant mechanism of seg- 
regation is a lateral sliding of the polymers with negligi- 
ble transverse displacements, then the linear density in 
Region I remains fixed at </>, defining the monomer den- 
sities (/)/ — 0, (/)// = 4>jii = '^^Zf ■ The different den- 
sities in regions I and II (HI) imply different blob sizes 
(Fig. mb)). Using Eq. ^ we obtain Fo = F{l,W,(f>l), 
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). For A > 1, 



f3F{l) 



and for < A < 1 



if 2i?|| -L<l<R\\ 
I] if < / < 2i?|| - L 



(4) 



(5) 



For A < 0, (3F{1) oc is identical to the free energy of 
two chains in an open tube [l^ , and the longitudinal con- 
finement has no effect. For A > 1, the total free energy 
is independent of I, therefore, there is no force driving 
segregation and the motion of the two chains is purely 
diff'usive For A < 1, F{1) is a monotonically in- 

creasing function of I, and this repulsive potential drives 
segregation (lit . 

The blob-structure-based free energy function (Eq. [5]), 
can be used to model the long-time behavior of the seg- 
regating polymers. The segregation time r is defined as 
the average time for chains moving from Z = to / = 0, 



which is equivalent to the mean first-passa ge t ime of the 
corresponding Fokker-Planck equation ©Mill. 



dP{l,t) 
dt 



-I3F{1) 



dl 



P{l,t) 



(6) 



where P(/, t) is the probability of the two chains overlap 
ping by / at time t, and D the diffusion constant given 
by the Einstein relation, D 
monomer friction coefficient 
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^ with ( being the 
3|. This description 
is valid if the segregation process is sufficiently slow such 
that at each stage, the motion of segments is controlled 
by equilibrium statistics [l^, [2^ . The global relax- 
ation time of the end-to-end distance of a confined poly- 
mer grows less rapidly than 
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As shown below, 
the segregation time in the regime A ~ 1, grows faster 
than A'^^, making the equilibration assumption reason- 
able. Since F{1) is a piecewise continuous function with 
different expressions in Regions I and II (III), the mean 
first-passage time can be written as the sum of the first- 
passage times of two subprocesses: t = ti + T2, where ri 
is for the subprocess of separation between I = and 
I — 2i?|| — L due to the entropy, and T2 is for separa- 
tion between I = 2i?|| — L and I = due to the competi- 
tion between the entropic force and longitudinal pressure 
2^ 2^- Introducing x = R\\ — I, we obtain, modulo con- 
stant factors: 
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where (31Fi{x) 



yX, j3T2{x) 



-^^'^yUydx (7) 
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In the limit of A « (L « 2i?||), the segregation time 
is governed by ti : 
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9f » 1 is 



In the second equation, the limit 
taken. Comparison with the time scale for pure diffu- 
sion from overlapping to separation r^iff ~ P^NR"^^ [B] , 
shows that the segregation time is much smaller than 
Tdiff. In the limit of A ~ 1, T2 is dominant. Since (3J^2{x) 
is a gradually decreasing function with changes of the or- 
der of (1 — A)^, a linear approximation can be used to 
estimate f3T2{x): 
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(9) 



Blob scaling arguments for free energy are only valid for 
(L — R\\) > W, which provides an upper limit: r < 
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the 



^^^^^^ ~Tdiff. Therefore, as L ^ K^^ , t -^Uiff, 
segregation dynamics approaches the unbiased-diffusion 
limit. Since ti, T2 are decreasing and increasing functions 
of L — , respectively, the segregation time can exhibit 
a minimum as a function of L for a fixed N and W. 

These theoretical predictions, were tested using Monte 
Carlo simulations based on the bond fluctuation model 
(BFM) Hi]. The BFM is a coarse-grained model 
of polymers in which chains live on a hypercubic lat- 
tice and fluctuations on scales smaller than the lattice 
constant are suppressed. The polymer is represented 
by a chain of effective monomers connected by bonds 
which are constructed to account for excluded-volume 
effects. An overlapping configuration of two chains is 
created by introducing a pseudo harmonic interaction: 
J2iLi ^(1^11 ~ R-2i)^, where R„m denotes the position 
vector of the nth monomer on the mth polymer and k 
is a parameter controlling the attractive strength. This 
interaction is turned off after the two chains are fully re- 
laxed in the overlapped configuration. The separation 
between two chains is measured by the horizontal {Xcc) 
and vertical (Ycc) projection of the centers of mass, and 
Xrr is related to I: 




if < ; < 2R\i - L 
if 2i?|| -L <l <R\\ 



(10) 



Simulations were performed over a wide range of pa- 
rameters: 80 < TV < 200, 30 < L < 140, keeping W = 10 
fixed. Each Monte Carlo trajectory spans a few hundred 
Rouse times 16|,|28|, and twenty independent trajectories 
are sampled for each set of parameters. We have checked 
the validity of the assumption in Eq. [3] and find that the 
difference of average extension between single chain and 
two chains is less than 6% [s^l- Figl^] illustrates the evo- 
lution of Xcc and Ycc with respect to Monte Carlo Steps 
(MCS). 

As shown in Figl^a), for A = 0.39,0.65,0.66, Xcc 
grows and fluctuates around L/2 and Ycc decreases to 
zero, a signature of segregation. The time scale for reach- 
ing a well-defined average is the same for Xcc and Ycc, 
which is a convincing argument for equating this mea- 
sured time to the calculated segregation time, r. Figl^j^b) 
shows that for A = 1.08, Xcc and Ycc do not grow but 
fluctuate between and [Xcc), and between and 
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(Fee), indicating a lack of segregation. 



In Fig. [2Ia),(b), Monte Carlo simulation results and 
numerical integrals of Eq. [7] are shown; both demonstrate 
non-monotonic behavior of r as a function of L. In Fig. 
[3Jc),(d), the segregation time r is scaled according to 
the relations in Eq. [8]and[9l the hnear fits support the 
theoretical predictions in the asymptotic regimes. 

Fig. m shows simulation data and the theoretical pre- 
dicted phase boundary separating the segregated and the 
mixed phase, demonstrating the accuracy of the theoret- 
ical predictions. Along each line, which corresponds to a 
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FIG. 2: (Color online) (a). Plot of 2X^^/L vs. MCS in the 
segregation region for (A = 0.39; iV = 100, L = 60), (A = 0.65; 
N = 120, L = 65), and (A = 0.66; iV = 200, L = 100). The 
inset shows Ycc vs. MCS. (b). Plot ol2Xcc/L and 2Ycc/L vs. 
MCS in the mixing region: (A = 1.08; = 100, L = 40). 



given chain length, the geometry with minimal segrega- 
tion time is marked (Fig. 3]) . A minimal extension of the 
blob picture to Zd ring polymers such as the chromosome 
of E. coli is to require ^biob £ ^ rather than W [s^l ■ Ap- 
plying this model to E.coli, using measured parameters 
EEallll: 6iob ^ 87nm, W 0.24nm, and L w 1.39/ito, 
locates E. coli in the segregation phase and close to the ge- 
ometrical condition of minimum segregation-time. While 
this observation may be fortuitous since chromosome 
strand is immensely more complicated than a linear poly- 
mer, it raises the interesting possibility that genome seg- 
regation times could have applied evolutionary selection 
pressure to genome lengths. 

In conclusion, we have investigated the dynamics of 
segregation of two self-avoiding chains confined to a nar- 
row rectangular box. A theoretical framework, based on 
the blob picture to capture the essence of the competition 
between excluded volume and confining effects, predicts 
a rich phenomenology of transitions between segregated 
and mixed states and optimal geometries that minimize 
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sions with Suckjoon Jun, Mike Hagan and Jane Kondev. 



FIG. 3: (Color online) (a) Plot of ri, T2 and r vs. L obtained 
from numerical integration of Eq. [7] with A'' = 80, W = 10. 
(b) Monte Carlo simulation results for the same set of parame- 
ters as (a), (c) Plot of r/iV^ vs. 1/(L-_R||) for 0.8 < A < 0.98. 
(d) Plot of t/N vs. (L - 7?||) for 0.03 < A < 0.2. In (c) and 
(d), the lines denote two fits for the asymptotic theoretical 
predictions (Eq. [HI [U . 
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FIG. 4: (Color online) Phase diagram of segregation and 
mixing. The x-axis parameterizes the monomer concentra- 
tion and the y-axis parameterizes the geometry. Simulation 
data are denoted by dots (segregation) and triangles (mix- 
ing). Lines illustrate data from various L: N — 80 (black), 
100 (red), 120 (blue), 150 (cyan), 180 (magenta), and 200 
(navy) from top to bottom. Stars mark aspect ratios with 
the minimal segregation times. The purple square denotes an 
estimate based on experimental data for E.coli. 



the segregation time. Monte Carlo simulations provide 
broad support for the theoretical predictions. Experi- 
ments in microfluidic devices should be able to provide 
direct tests of the predictions, and elucidate the role of 
entropy in driving segregation of biopolymers. The sim- 
ulation can be extended to study more realistic models 
based on actual chromosomes structure of bacterial. 
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